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ABSTRACT 

In this paper we present the results of a survey for radio transients using data obtained 
from the Very Large Array archive. We have reduced, using a pipeline procedure, 5037 
observations of the most common pointings - i.e. the calibrator fields. These fields 
typically contain a relatively bright point source and are used to calibrate 'target' 
observations: they are therefore rarely imaged themselves. The observations used span 
a time range ~ 1984 - 2008 and consist of eight different pointings, three different 
frequencies (8.4, 4.8 and 1.4 GHz) and have a total observing time of 435 hours. We 
have searched for transient and variable radio sources within these observations using 
components from the prototype LOFAR transient detection system. In this paper 
C*") . we present the methodology for reducing large volumes of Very Large Array data; 

and we also present a brief overview of the prototype LOFAR transient detection 
algorithms. No radio transients were detected in this survey, therefore we place an 
upper limit on the snapshot rate of GHz frequency transients >8.0 mJy to p ^ 0.032 
deg -2 that have typical timescales 4.3 to 45.3 days. We compare and contrast our 
upper limit with the snapshot rates - derived from either detections or non-detections 
r> ■ of transient and variable radio sources - reported in the literature. When compared 

with the current Log N - Log S distribution formed from previous surveys, we show 
that our upper limit is consistent with the observed population. Current and future 
radio transient surveys will hopefully further constrain these statistics, and potentially 
discover dominant transient source populations. In this paper we also briefly explore 
the current transient commissioning observations with LOFAR, and the impact they 
will make on the field. 
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1 INTRODUCTION 

The scope and depth of transient radio science is vast: by 
utilising the time domain we can gain unique insight into 
such objects as neutron stars and white dwarfs in binary 
systems, relativistic accretion and consequent jet launch 
around black holes, distant gamma-ray burst afterglows, su- 
pernovae and active galactic nuclei (AGN), to name a few. 



The distances to these objects, as well as the timescale for 
transient behaviour, varies dramatically. For example, giant 
kilo-Jansky micro-second radio pulses have be en observed 
from the relatively nearby Crab Pulsar (e.g. see Bhat et"ail 
2008). In contrast, month timescale (and longer) variations 
are commonly observed in the radio emission produced by 
powerful jets driven b y accretion onto super-massive black 
holes in distant A GN (|de Vries et al.ll2004l . iBell et alJliuTol . 
IJones et al.ll2010T ) . Through studying the transient and vari- 
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able nature of these exotic and energetic objects, we obtain 
an unprecedented laboratory to probe extreme physics. 

Despite the scientific potential, the transient and 
time variable sky is a relatively unexplored region 
of parameter space. Historically radio transient de- 
tections have been sparse due to an inefficient sur- 
vey figure of merit to adequately sample a large 
amount of sky to sufficient sensitivity a n d time reso- 
lution l|Cordes. Lazio. fc McLaughlin! |2004| ; iHessels et al.l 
2009). Some detections of transients have therefore 
been made serendipitou sl y (for examples see iDavies et al] 
19761 . IZhao et al. | Il992l. Ivan den Oord fc de Bruvnl 1 1994 



Bower et al .1120031 and 



Lenc et alj|2008). This limitation will 



soon be relieved by the next generation of wide field tele- 
scopes and their respective dedicated transient surveys. 

A variety of new wide-field facilities will soon be avail- 
able to sample the transient sky. I n the optical ban d the 
Palomar Transient Factory (PTF; iRau et al] 120091 ) and 
Panoramic Sur vey Telescope and R apid Response System 
(Pan-STARRS; lHodapp etZll2004l ) will survey the sky for 
transie nts. In the radio band the Allen Telescope Array 
(ATA; I Welch et all 120091'). the M urchison Wide Field Ar- 
ray (MWA; lLonsdale et al] 20091) an d the Low Frequency 



Array (LOFAR; Fender et al. 20081 ') will soon begin or 



have already commenced operations. Other wide field radio 
pathfinders such a s the Karoo Array Telescope (MeerKAT; 
iBooth et al] I2009T ) and the Australian Square-K ilometer- 
Array Pathfinder (ASKAP: [Johnston et al]|200Sl ) are also 
being developed on the road to the Square Kilometer Array 
(SKA). Transient studies are a key science goal for all of 
these facilities. 

A common method to detect radio transients is 
through multi-wavelength triggered observations from, for 
example, all sky monitors on X-ray observatories. These 
have produced radio counterparts to gamma-ray burst 
(GRB) afterglows and black hole X-ray binary outbursts 



(for examples see Frail et all 1 19971 . iGaensler et al] 120051 . 



lEck. Cowan, fc Branch! l2002h . This method relies on hav- 
ing a detectable high frequency counterpart, which may be 
absent (or difficult to detect) for sources such a s X-ray dim 
isolated neutron stars (XDINs; Ofek et al. 20091) and orphan 
gamma-ray burst afterglows (| Frail et al] 19971 ) , demonstrat- 
ing the need for dedicated radio transient programs. 

Despite the historical challenges, dedicated or commen- 
sal transients surveys have produced a number of interesting 
results. The Galactic centre (GC) has been the area for some 
intense observing campaigns; these studies have so far de- 
tected a number of radio transients, th e most recent being 
GCR T J1742-3001 (iHvman et al] 120091) and GCRT J1 746- 
2757 <|Hvman et aTll2002l ) - also s eelBower et al] (120051). In 
the high time resolution domain, iMcLaughlin et al] (|2006l ) 
discovered short duration transient radio bursts from neu- 
tron stars, i.e. RRATs (Rapidly Rotating Transients) - as 
well as many detections of new pulsars (for example, from 
the Parkes multi-beam survey). 

Nine bursts - named the WJN transients - in ex- 
cess of 1 Jy have been discovered using drift scan obser- 
vations with the Was eda Nasu Pulsar Obser vatory at 1.4 
GHz (summarised in iMatsumura et aT 2 009 | . but also see 



Kunivoshi et a JlgOOirNiinuma et al.ll2007l.lKida et al.ll20o"sl. 



Niinuma et al . 2009] for further details). These are some of 
the brightest transients reported in the literature and so far 



remain unexplained. Recently [Croft et al] (|2010h published 
results from the ATA Twenty Centimetre Survey (ATATS): 
no transients were detected and an upper limit on the snap- 
shot rate of events was given. Subsequently the Pi GHz Sky 
Survey (PiGSS) surveyed the sky with the ATA at 3.1 GHz, 
providin g the deepest static source catalogue to date above 
1.4 GHz jBower et al]|2010l) . No transients were reported in 
this survey and an upper limit on snapshot rate was placed. 

Searching for highly variable known radio sources can 
also be a useful diag nostic in examining the dynamic ra- 
dio sky. For example, ICarilli et al] (|2003l ) found a number 
of highly variable (A5* ^ ±50%) radio sources in a small 
number of repeated observ ations of the Lockman Hole at 
1.4 GHz. iFrail et al] (|2003l ) also found four highly variable 
radio transient sources from follow-up observations of GRBs 
at 5 and 8.5 GHz: the rates of these event s are consistent 
with those reported in Carilli et al] 1200 ^). Reporting 39 
variable radio sources, iBecker et al] (|20 lOf) recently charac- 
terised the surface density of variables in the direction of the 
Galactic plane at 4.8 GHz. Most of the variable sources de- 
tected had no known multi-wavelength counterparts. This 
is an important study as the rates of transient and vari- 
able sources may differ in the direction of Galactic plane 
when compared with an extra-Galactic pointing. In particu- 
lar, large numbe rs of flare stars are known to produce bright 
coherent bursts l) White et al.lll989 | ). and could do minate de- 
tections at low frequencies (jBastian et al"1 [l988). A deeper 
discussion on the difference between transient and variable 
processes will follow later. 

Radio telescope archives potentially contain many hours 
of data which currently remains unsearched for radio tran- 
sients. An arc hival study comparin g the NVSS (NRAO VLA 
Sky Survey; ICondon et ail 1 19981 ) and FIRST (Faint Im- 
ages of the Radio Sky at Twenty-cm; IBecker et aD Il995l : 
IWhite et al] Il997n catalogues was conducted by Levin- 
son et al. (2002), with a follow-up study by Gal- Yam et 
al.(20 06): a numbe r of r adio transient sources were iden- 
tified. Bower et al] (|2007f ) analysed 944 epochs of archival 
VLA data at 4.8 and 8.4 GHz spanning a period of 22 
years. In this survey ten radio transients were reported, 
with the host galaxies possibly identified for four out of the 
ten sourc es, and the hosts and p rogenitors of the other six 
unknown. [Bannister et ail l|201ll ) recently published results 
from a search for transient and variable sources in the Mo- 
longlo Observatory Synthesis Telescope (MOST) archive at 
843 MHz. 15 transient and 53 highly variable sources were 
detected over a 22 year period. iBannister et al] (|201lh use 
these detections to place limits on the rates of transient and 
variable sources. iBower fc Sau 3 (|2010D have published fur- 
ther archival work examining observations of the calibrator 
3C286 at 1.4 GHz. A total of 1852 epochs are examined 
spanning a time range 23 years: no radio transients were 
reported. 

IBower et al] l|2007l ) with its high yield of transients pro- 
vided the motivation for this work. The aim of this paper is 
to push archival radio transient studies further, within the 
framework of testing and refining the transient detection al- 
gorithms that will operate on the LOFAR radio telescope 
|Fender et~aT]|2008t ). We present the findings of an 8.4, 4.8 
and 1.4 GHz study of the repeatedly observed flux and phase 
calibrator fields found in the VLA archive. We explore from 
some of the publications discussed above the reported snap- 
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Table 1. Number of images reduced and searched with respect to observing frequency. Note, 4.8 and 8.4 GHz were the primary frequencies 
of interest, 1.4 GHz images were only produced for one field. The calibrators are referred to by their J2000 epoch name. < 5T nex t > is the 
average change in time between sequential observations (including observation on the same day - see section 2.1). <t> is the average 
integration time spent on the calibrator. 



Field 


R.A. 


Dec. 


t 


b 


# Obs. 


# Obs. 


# Obs. 


< &T nex t > 


<T> 




(J2000) 


(J2000) 






(8.4 GHz) 


(4.8 GHz) 


(1.4 GHz) 


(days) 


(mins) 


1800+784 


18 h 00 m 45 3 .7 


+78°28 m 04 s .0 


110.0 


+29.1 


992 


908 


151 


4.3 


5.6 


0508+845 


05^08*" 42 s .4 


+84°32 m 04 s .5 


128.4 


+24.7 


205 


413 




13.6 


5.1 


1927+739 


19' l 27 m 48 a .5 


+73°58 m 01 s .6 


105.6 


+23.5 


171 


19 




45.3 


9.0 


1549+506 


15 fe 49 m 17 s .5 


+50°38 m 05 s .8 


80.2 


+49.1 


480 


558 




8.3 


5.0 


0555+398 


05 h 55 m 30 s .8 


+39°48 m 49 s .2 


171.6 


+7.2 


123 


190 




27.3 


3.5 


2355+498 


23' l 55 m 09 a .5 


+49°50 m 08 s .3 


113.7 


-12.0 


183 


99 




29.3 


7.6/15.3 a 


3C48 


0l' l 37 m 41 a .3 


+33°09 m 35M 


25.1 


+33.4 




545 




16.2 


4.6 


Total 










2154 


2732 


151 


1.9/3.0 6 


5.2 



a 7.6 mins at 8.4 GHz and 15.3 mins at 4.8 GHz. 

b 1.9 days including observations on the same day / 3.0 days disregarding them. 



shot rates of either detections of radio transients, or upper 
limits based on non-detections: we place our own upper limit 
on these values and discuss the implications. 



2 VLA DATA 

2.1 VLA calibrator fields 

We have searched part of the VLA archive for transients. In 
order to optimise our chances of success we searched for the 
most repeatedly observed fields in the VLA archive. The 
flux and phase calibrator fields were chosen as the back- 
bone of this new transient study. These calibrators are ob- 
served routinely and are a standard and necessary calibra- 
tion technique in radio interferometry. The calibrator fields 
usually contain a relatively bright compact object, typi- 
cally a quasar. The calibrators were selected to fulfil one 
or more of the following criteria: (A) they were observed 
frequently; (B) they should be unresolved on the longest A- 
configuration VLA baseline; (C) they had a relatively large 
integration time per observation. The chosen VLA fields are 
summarised in Table 1. Note that we initially focused our 
efforts on the flux calibrators, specifically 3C48; however, 
with the nature of the bright source in the field and the typ- 
ically short integration time (~ 1 minute) the images often 
suffered from artifacts. Therefore we quickly switched our 
attention to the phase calibrators - which proved, due to 
longer integration times (~ 5 mins) to have better image 
fidelity. 

A total of 5037 flux and phase calibrator images have 
been reduced and searched, with a total observing time of 
435 hours. The average integration time spent on all sources 
was < t > = 5.2 mins. A full statistical description of the 
measured noise per image is given in section 5. 

The mean separation between observations regardless of 
pointing and frequency was < ST next i > = 1.9 days. When 
calculating this average we have included observations that 
occurred on the same day. Note that it is quite common 
for an observation of the same source to be performed at 
a number of different frequencies (i.e. sequentially within 
the same observation). This obviously produces a bias that 



reduces the average time between observations. When pro- 
ducing the images for this survey we only logged the date of 
the observations, not the exact start and stop time. Extract- 
ing the start and stop time was not easily executed within 
the imaging pipeline framework. We therefore add in an ar- 
bitrary time delay of 4 hours in calculating the averages for 
observations that have the same date. Ignoring observations 
on the same date yields an average < 5T next 2 > = 3.0 days 
(regardless of pointing and frequency). 

See the top panel of Figure[T]for a histogram of the time 
differences between sequential observations for all point- 
ings and frequencies. The bottom panel of Figure [1] shows a 
(summed) histogram of the time differences between obser- 
vations per pointing. Table 1 summarises the average time 
difference between observations and the average integration 
time per pointing. The lowest average time difference be- 
tween observations achieved for one single pointing was 4.3 
days (1803+784); the highest was 45.3 days (1927+739). As 
we have sampled a number of different fields at different ca- 
dences we state the range 4.3 to 45.3 days as the timescale of 
transient behaviour that we would be sensitive to. In quoting 
these numbers we do make the assumption that each image 
has an equal chance of making a detection - i.e that there is 
an isotropic distribution of radio transients and there is no 
frequency dependence (between 1.4 and 8.4 GHz) for detec- 
tion. 



2.2 VLA data analysis 

The archival data have be en reduced using the ParselTongue 
(see iKettenis et al.l 120061 ) Python interface to the Astro- 
nomical Image Processing System (AIPS; iGreise 3 l2003h . 
A scripted procedure was written to perform the following 
tasks: 

• The data were loaded into AIPS. 

• The antenna table was searched for antennas that were 
'out' or designated 'EVLA' (only relevant 2006 onwards) 
during the observations. These antennas were then flagged 
from being used in subsequent calibration and imaging 
tasks. 

• The maximum baseline length in the observation was 
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Figure 1. Top Panel: Histogram showing the time difference be- 
tween observations - for all pointings and frequencies within this 
survey. Bottom Panel: A (summed) histogram of the time differ- 
ences between observations per pointing including all frequencies. 
For a description of the average cadence per pointing see table 1. 



ascertained; this was then utilised to optimise cell and image 
sizes respectively in the subsequent imaging steps. 

• The automated flagging procedures were applied to flag 
unwanted and erroneous visibilities. 

• The observation log was searched for any known flux 
or phase calibratoiQ used at the VLA within an entire ob- 
servation. 

• For any identified calibrator, standard calibration was 
performed. 

• The calibrator was subsequently imaged and decon- 
volved, followed by three iterations of phase self-calibration 
only, and then one iteration of amplitude and phase self 
calibration. The time interval for phase calibration was set 
appropriately with respect to the integration time on source. 

• The calibrator was boxed off, modelled and removed 
from the image using UV component subtraction. 

• Finally the su btracted image was lightly CLEANed 
with 150 iterations l|Hogbomlll974T ). 

The script was designed to be run on large volumes of 
data without interruption. Python exception handling was 



1 Found at http://www.aoc.nrao.edu/~gtaylor/csource.html 



used to catch potential errors and remove bad data from 
further processing. An example image of 3C48 produced by 
the pipeline before and after source subtraction is shown in 
Figure 2. Both images show contours and Grey scale to give 
the reader an intuitive feel of the image quality. The source 
to the North of 3C48 is persistent with a flux ~ 30 mjy and 
will be discussed further in section 4. 



3 TRANSIENT SEARCH 

The images produced by the imaging pipeline were then pro- 
cessed thr ough the prototy pe LOFAR transient detection al- 
gorithms (ISwin bank 20071) ■ By prototype, we refer to a well 
tested subset of algorithms (taken from LOFAR) that were 
put together in a 'pipeline' to perform the source extraction 
and databasing for this transient search only. It should be 
noted that a broader set of algorithms with greater function- 
ality will define the final LOFAR transient detection system. 

For each image a background RMS map was calculated 
over the entire image. For any island of pixels above 8cr i.e. 
eight times the noise measured from the RMS map, source 
extraction was performed by fitting elliptical Gaussians. For 
further details on the LOFAR source extraction algorithms 
see Spreeuw (2010). In the case of LOFAR images, a sub- 
section of the entire image will be used to calculate a 'local' 
RMS map. Considering the increased field of view of LO- 
FAR, coupled with greater source counts, it is not sufficient 
to define a global RMS map over the entire image. 

A conservative threshold of 8a was chosen after initial 
tests and pipeline refinements indicated that a badly cali- 
brated and reduced image - considering the typical strength 
of the calibrator flux - could produce a large number of arti- 
facts and thus false source detections. All images that were 
processed by the imaging pipeline were processed through 
the prototype LOFAR transient detection algorithms - er- 
rors included. The rationale for this was to explore the effects 
of badly calibrated images with respect to source extraction 
and transient detection. LOFAR will incorporate a false de- 
tection rate (FDR) algorithm in the source extraction sys- 
tem, where the global detection threshold for source extrac- 
tion is set to minimise the number of false p ositives and is 
gover ned by the individual image statistics (jHopkins et al.l 
120021 ). 

After source extraction a MonetDB (Boncz 2002) 
database was then populated with the measured properties 
and associated data of the extracted sources. The source 
properties include: position and associated errors, all Stokes 
parameters of peak and integrated flux including the Gaus- 
sian fitting parameters. The associated data included, for 
example, time of observation, operating frequency and beam 
properties. The database was searched for either unique 
sources that had no previously known counterparts (in pre- 
vious images), or a known source whose flux had varied 
by a significant amou nt. We adopt the same metric used 
in ICarilli et al.l l|2003l ) to define significant variability as 
AS > ±50%. 

For further details on the transient database algo- 
rithms, including a mathematical description of the vari- 
ability measure used in this survey see Scheers (2010). The 
database will also cross-reference the transient parameters 
with those in the WENSS (Westerbork Northern Sky Sur- 
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Figure 2. An example image of 3C48 at 4.8 GHz produced by the pipeline (observation date 1984-10-05). The left panel shows a 
CLEANcd image before source subtraction. The right shows the same epoch after source subtraction, a cross denotes the original 
position of 3C48. The wedge on the right of each image shows the intensity of the Grey scale in Janskys. The integration time was ~4 
mins, yielding an RMS of 0.6 mjy. Contour levels are -3, 3, 8, 20, 40, 100, 500, 2000, 5000, 8000 X RMS for both images. 



vey; iRengelink et~afl 1 19971). NVSS and VL SS (VLA Low- 
frequency Sky Survey; I C ohen et al.l [2007) catalogues to 
search for counterparts; for this survey we typically relied 
on NVSS. 

For the majority of fields in this survey, after the cal- 
ibrator source was subtracted, the fields were left almost 
devoid of sources: thus the transient search was relatively 
trivial. A few ~mjy radio sources were present in some of 
the fields and only some of the time, due to changing sen- 
sitivity. Therefore, we concentrated our efforts on locating 
unique sources, rather than characterising the variability of 
known sources. Once the final list of candidate transient 
sources had been produced, light curves were automatically 
generated and the images of interest were checked for cali- 
bration errors and image fidelity. 



4 RESULTS 

A total of 5037 images at various pointings and frequencies 
have been searched at a detection level of 8a. Nine candidate 
transients that were detected in images with adequate image 
fidelity were scrutinised further. Four of these candidates 
were found to lie consistently on the dirty beam. After re- 
reduction they were shown to be calibration errors which 
had been cleaned to a point source. 

One candidate was detected 76 times at 4.8 and 8.4 GHz 
and was found to be significantly variable. After careful con- 
sideration and review of the literature it was concluded that 
this source was created by a bug in the VLA recording sys- 
tem, whereby the pointing of telescopes was changed with- 
out updating the header information; this error had affected 
a p reviously reporte d transient VLA J172059. 9+385226. 6 
(see lOfek et ai]|2009l for further details). This error was not 
detected in any of the other calibrator images taken around 
the same time. 

Three of the candidates were found to be associated 



with known weaker radio sources. These sources were de- 
tected in the deepest observations (~ 30 mins) of the respec- 
tive fields. As observations at this depth were very sparse, 
these radio sources were considered - by the database al- 
gorithms - as transient. These candidates were quickly re- 
moved when cross referenced with the radio catalogues. The 
last transient candidate although surviving some of the re- 
reduction tests was discarded after it was re-calibrated and 
the significance level dropped below adequate levels. 

The 3C48 field did contain a persistent source at a = 
01 h 37 m 44 s 2 and s _ + 33°ii™26 s .5 (J2000) with S v ~ 30 
mjy. This source could not be identified in the NVSS cata- 
logue, due to insufficient resolution to separate it from 3C48. 
The FIRST survey did not cover the position needed to cat- 
alogue this source. The source is however pre viously ident i- 
fied in high dynamic range studies of 3C48 fsee lBrigdll995h . 
This source was searched for significant variability but none 
was found. 



5 SURFACE DENSITY UPPER LIMIT 

As we have detected no radio transients with this survey, 
we use the area surveyed per observation, coupled with the 
typical sensitivity to constrain the snapshot rate of transient 
events. To calculate the 2a upper limit of the snapshot rate 
of transients from our survey we assume a Poisson distribu- 
tion; for zero detections (n=0) we use: 

P(n) = e' pN (1) 

where p is the snapshot rate of transients; and the 2a con- 
fidence interval is defined as P(n) = 0.05 at the 95% confi- 
dence level. N is the sum of the number of images multiplied 
by the field of view (fi) at that given frequency i.e. 

N = (n 8 . 4 x N SA ) + (fi 4 .s x N 4 .s) + x N 1A ) (2) 
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Figure 3. Top panel: Histogram showing the measured noise - 
calculated in the same region - for all images. < at > gives the 
theoretical noise calculated from the average integration time for 
all observations; < 10 X at > gives this limit multiplied by 10. 
Bottom panel: Histogram showing the measured noise divided 
by the theoretical per observation - which accounts for different 
integration times. 



Note, we only consider a search area within the half-power 
radius per image. Evaluating equations 1 and 2 yields a snap- 
shot rate of p ^ 0.032 deg -2 . 

To evaluate the flux density limit that we are sensitive 
to when searching for transients, we must statistically con- 
sider the noise measured in all the images. Figure 3 presents 
this information in two different ways. Firstly, on the top 
panel we show a histogram of the measured noise a m in all 
the images; included on the plot is an indicator of the the- 
oretical noise in a 5 minute observation (0.16 mJy at 4.8 
GHz), and also ten times this value. Secondly - in the bot- 
tom panel - we show a histogram of the measured noise di- 
vided by the theoretical noise a t (i.e. a m /a t ) for all images: 
thereby taking into consideration that not all observations 
were ~ 5 minutes in duration. 

It can be seen that the bulk of the images achieved an 
image noise less than Wat- Possible deviations away from 
the theoretical noise could be attributed to, for example, un- 
removed RFI, bad calibration solutions, effects of a bright 
source in the field and, in general, settings and assumptions 
within the imaging pipeline that do not lend themselves to 
a given observation. Taking the median value of all the mea- 



sured noises we find a me dian ~ 1 mJy (or 6.25<at >)• Using 
8a m edian as the global detection threshold for the entire sur- 
vey, we find that we would be sensitive to transients >8 mJy, 
with typical timescales 4.3 to 45.3 days (see section 2.1). 

In calculating this upper limit we have included all im- 
ages reduced by the pipeline; however, not all images were 
reduced successfully. Although some images contained arti- 
facts, they rarely contaminated the entire image, thus some 
area could still be searched effectively. If a unique transient 
point source was detected in an image with errors, the image 
was re-reduced by hand and checked for reproducibility (see 
section 4). 

In Figure 4 we compare the limit imposed on the snap- 
shot rate of sources from this study, with those found in the 
literat ure - this figur e is derived f rom F ig. 9 of lBower et al.l 
|2007l ) and Fig. 20 of ICroft et all (|2010l V We do not include 
a typical timescale of the transient duration in the snapshot 
rate calculation as it is not well constrained. This informa- 
tion is summarised and referenced in Tableland described 
further in the following paragraphs. 

The iBower et al.l (|2007f ) survey reported the snapshot 
rate of transients to be p=1.5±0.4 deg -2 (labelled 'B2007 T 1 
Week' in Figure 4) from eight detections, with characteristic 
timescale 20 minutes < t c h ar < 7 days, above a flux density 
370 pJy (with typical image noise ~ 50 pJy at the pointing 
centre). Two transients were detected in the 2 month av- 
eraged images above a flux density of 200 /iJy, giving a 2a 
limit on the snapshot rate p ~2 deg -2 (labelled 'B2007 T 2 
Month' in Figure 4). No transients were detected in the year 
long averages above 90 pJy: limiting the 2a snapshot rate 
to p <6 deg" 2 (labelled 'B2007 T 1 Year' in Figure 4). 

The PiGSS-I survey using the ATA at 3.1 GHz sets 
an upper limit on the snapshot rate of transients to p < 1 
deg -2 at 1 mJy, and p < 0.3 deg -2 at 10 mJy (labelled 
'B2010 T ' V ' in Figure 4), with characteristic timescale one 
month l|Bower et al.ll20Toh . A recent study bv lBower fc Saull 
i|2010T) - using archived VLA observations of the flux cali- 
brator 3C286 at 1.4 GHz - set an upper limit of p < 3 x 10 -3 
deg -2 at 70 mJy, and p < 9 x 10~ 4 d e g" 2 a t 3 Jy (labelled 
'B2011 T '). The work of lBower fc Saull (|2010T l is very compa- 
rable to the work presented in this paper, however, in this 
study by using predominantly the phase calibrator fields, we 
slightly push the mean sens itivity down. 

The lCarilli et all (|2003l ) study set an upper limit on the 
rate of highly variable radio sources ^ 100 pJy to < 18 deg -2 
with characteristic timescales of 19 days and 17 months (la- 
belled 'C2003 v ' in Fi gur e 4). Note that the se were detections 
of variable radio sources. iFrail et al.l (|2003h derived a compa- 



rable quantity to lCarilli et all |2003l) of p ~ 5.8 deg" 2 with 
four highly variable sources with characteristic timescale ~ 1 
day, above a flux density 250 pJy (labelled ^2003^ ' in Fig- 
ure 4). Similar to these surveys bu t in the direction of the 
Galactic plane iBecker et al.l Q2010) found 39 variable radio 
sources between a flux density range 1 to 100 mJy, varying 
on timescales of years: they derived p ~ 1.6 Galactic sources 
deg" 2 (labelled 'Bk2010 v ' in Figure 4). 

In the context of this survey, the difference between a 
variable and a transient - from a purely observational sense 
- is a matter of detectability: transients sit, most of the time, 
below the detection capabilities of the instrument; while 
variables sit above or close to it. However, the underlying as- 
tronomical processes associated with transient and variable 
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Becker et al. (2010) 
This work 
Bower et al. (2007) 
Bower et al. (2010) 
Bowers Saul (2010) 
Carilli at al. (2003) 
- Croft etal. (2010) 
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Figure 4. Snapshot rate (deg -2 ) against flux density (Jy) of detections of transients (labelled 'T'), detections of variable sources (labelled 
'V') and upper limits based on non-detections (labelled with downward arrows). The thick black line denotes detections; the thin line 
denotes upper limits. The Bower et al. (2007) 1 week, year and two month limits are indicated as B200 7 r with the appropriat e time-scale. 
Ban2011 v and Ban2011 T indicates the separate rates derived for variables and transients reported in lBannister et all l|201lh . Bell2011 T 
indicates the 2rr upper limits derived from this study. 'LOFAR F ' indicates the theoretical constraint that LOFAR could provide with 
zero detections from 10 epochs of 12 hour observations, each of 25 deg 2 fields (using 18 core and 18 remote stations at 150 MHz); 
'LOFAR c ' indicates the current commissioning capabilities at 150 MHz. 'LOFAR F 30 MHz' shows the rate calculated for a 30 MHz field 
of view (assuming the final theoretical noise is reached with 18 core and 18 remote stations). We note that this plot does not contain 
any information on characteristic time duration and recurrence of transient behaviour as both are currently poorly constrained. 



sources may differ, and should require a different treatment 
when considering the rates of events. For example, we might 
expect the rates of variable sources to differ from 'one off' 
explosive transients such as GRB afterglows - which have a 
finite lifetime and will be undetectable beforehand. For fu- 
ture surveys, a spectrum of transient and variable behaviour 
will be observed depending on the cadence and sensitivity. 
The boundaries between the definitions will become more 
blurred as the cadence and sensitivity is increased. 

The WJN transients summarised in iMatsumura et all 

(2009) range in flux density from 1 to 4.3 Jy with character- 
stic timescale ~ 1 day, yielding a snapshot rate p ~3xl0 -3 
deg" 2 (labe l led 'M at2009 T ' in Figure 4). In comparison the 
ICroft et alj (|2010h survey set a 2<r upper limit on the snap- 



shot rate of events > 40 mjy to be p <0.004 deg -2 ; by 
comparing their source fluxes with those in the NVSS cat- 
alogues the characteristic timescale is ~ 15 years (labelled 
'Croft2010 T 'in Fi gure 4) . The most st ringent limit placed on 
the snapshot rate of sources is set bv lGal-Yam et all (|2006l ) 
to be p < 1.5 x 10 -3 deg -2 for flux densities >6 mjy (labelled 
'GY2006 T ' in Figure 4). Note, the FIRST survey has im- 
proved angular resolution (5") when compared with NVSS 
(45"), therefore correct source association effects transient 
identification. We do not state a characteristic time scale 
for the FIRST-NVSS comparison as both individual surveys 
took a number of years; specific timescales can only be con- 
sidered on a source by source basis. 

iBannister et al.l (|201 lh set a limit on the snapshot rate 
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Table 2. Summary of snapshot rates reported in the literature. The results are separated out according to upper limits based on non- 
detections (top); transient detections (middle); and detections of highly variable radio sources (bottom). The flux min column designates 
the detection threshold of the observations report ed in the literature or the minimum flux of detections (indicated as such); the maximum 
flux is only indicated for transient detections. The lBower et all ll2007fl results have been stated three times depending on the characteristic 
time s cale sampled. We do not g ive the number of epochs for the WJN transients as it is not stated clearly in the literature. iBower et al.l 
l|2010l) andU ower fc Saull (l2010f) state two different rates depending on flux density, we quote these separately as (A) and (B). 
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0.0738 
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Bannister et al. (2011) 


14000(5er) 


6.5xl0 6 


1.3 x 10~ 2 
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3011 6 


Bower et al. (2007) 


370 


7042 


1.5±0.4 


20 mins - 7 days 


4.8 and 8.4 


944 


Bower et al. (2007) 


200 
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2 


2 months 


4.8 and 8.4 


96 


WJN/Matsumura et al. (2009) 


lxlO 6 


4.3xl0 6 


3xl0~ 3 


~ 1 day 


1.4 




Bannister et al. (2011) 


> 14000 




0.268 


days - years 


0.843 


3011 6 


Carilli et al. (2003) 


>100 




<18 


19 days and 17 months 


1.4 


5 


Becker et al. (2010) 


>100 




1.6 


~ 15 years 


4.8 


3 b 


Frail et al. (2003) 


>250 




5.8 


~ 1 day 


5 and 8.5 





a Different noise values were found in each survey map thus global threshold taken above 6 mJy 
b Combined mosaic 



of transient sources (calculated from detections) at 0.848 
GHz to be p < 1.3 x 10 -2 deg -2 above 14 mJy at a variety 
of timescales (labelled 'Ban2011 T ' in Figure 4). For variable 
radio sources a rate of p < 0.268 deg -2 is expected between 
a flux density 14 to 100 mJy (labelled 'Ban2011 v ' in Figure 
4). 

The upper limit derived fro m this study i s con sistent 
with the detections reported by I Bower et al.l (|2007| ) - we 
might have expected possibly one detection at our flux den- 
sity thresholds, assuming that t he transient popula tion sam- 
pled is isotropically distributed. IBower et all (120071 ) did note 
that an overdensity of galaxies was found within their field. 
The rate derive d from this work i s also broadly consis - 
tent with that oflBower et al ] ||201(t).lBannister et al.l (120111 ) 



and IBower fc Saull (l201oT yif the 



Bow er et all (|2007f ) and 



iBannister et al. ( 20 111 ) detections are of a similar nature, 



then some of the surveys that report upper limits sit very 
close to the 'real' Log N - Log S. This is the best benchmark 
to date to predict the parameter space a given survey should 
probe to find transients. However, measurements such as fre- 
quency dependence and characteristic timescale of transient 
behaviour still need to be constrained. 



6 PREDICTIONS FOR LOFAR 

Commissioning observations are currently underway with 
LOFAR that are probing the parameter space described in 
this paper. We include in Figure 4 a currently theoretical up- 
per limit of p < 0.012 deg -2 based on zero detections from 



10 epochs of 12 hour observations, each of 25 deg 2 fields 
at 150 MHz (labelled LOFAR c ). Early observations with 
LOFAR around August - September 2010 yield a typical 
RMS of 15 mJy (75 mJy for a 5<r detection which is plot- 
ted) with a bandwidth of 50 MHz spread over 256 sub-bands 
(16 channels per subband). However as more baselines have 
come on-line, and next-generation data reduction strategies 
have been implemented, improvements upon this value have 
been made. A more realistic final theoretical noise of 0.36 
mJy based on 18 core and 18 remote stations is indicated in 
Figure 4 (labelled LOFAR F ). We also include a theoretical 
prediction based on 10 epochs of 419 deg 2 fields at 30 MHz 
(labelled 'LOFAR F 30 MHz' ). This sets an upper limit on 
the snapshot rate of transients to p < 0.0018 deg~ 2 above a 
detection threshold 30 mJy (6 mJy RMS noise). 

Recent work by lLazio et al.l |2010l ) using the Long 
Wavelength Demonstrator Array (LWDA) - a 16 dipole 
phased array with all-sky imaging capabilities - have per- 
formed an all-sky blind transient search. A total of 106 hours 
of data was searched for radio transients at 73.8 MHz - the 
largest survey yet (in imaging mode) at low frequencies. 
With no detections of radio transients outside of the solar 
system above a flux density 500 Jy, an upper limit of 10~ 2 
yr -1 deg -2 is placed on the rate of events. With a typical 
integration time of 5 minutes, this converts to p < 9.5xl0 -8 
deg -2 ; we plot this limit in Figure 4 (labelled 'L az2010') as- 
sumin g a 5cr (2500 Jy) detection is needed. The lLazio et al.l 
l|20ld ) results tell us that extremely bright radio transients 
with characteristic timescale ~ 5 minutes are very rare. Ob- 
servations with LOFAR at 30 MHz (see Figure 4) would be 
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complementary to the lLazio et all (|201Ch survey. Approxi- 
mately 20 tiled pointings could offer the same solid angle 
coverage as the LWDA i.e. the whole sky, with increased 
sensitivity. Pushing into this parameter space on a logarith- 
mically spaced range of timescales is a goal for LOFAR, 
as well as an all-sky monitoring functionality to catch the 
brightest and rarest exotica. 

We can see from Figure 4 that if LOFAR observations 
were separated ~ weekly, we wo uld be able - via sa mpling 
similar parameter space to the iBower et all (|2007f ) detec- 
tions - to test the differences between a GHz and MHz pop- 
ulation of radio transients. If the emission mechanism for 
the GHz population is predominantly via the synchrotron 
process, then many sources will be initially optically thick 
within the LOFAR band. The rise time for a distant, lumi- 
nous, population of radio transients - such as GRB after- 
glows - could be months or years; with lower peak fluxes. 
Therefore a steep spectrum population of coherent emitters 
might dominate detections in the LOFAR band. This coher- 
ent population will not be limited by the brightness temper- 
ature limit of the synchrotron sources and they could also 
have more erratic cadences - i.e. switch on and off - which 
should in turn affect the snapshot rate of events. 

Over ten observations of the same field, with approxi- 
mately a weekly cadence have already been obtained with 
LOFAR at 150 MHz. The data reduction is in progress and a 
concise transient search will follow soon. We therefore hope 
to test the hypothesis above shortly and push further using 
future experiments with LOFAR. 



7 CONCLUSION 

In this paper we have presented the results of an archival 
VLA study. We have calibrated, imaged and searched 5037 
images of the calibrator fields totalling 435 hours of observ- 
ing time. We have presented the methodology for reduc- 
ing VLA data in a pipeline procedure. We have also ex- 
plored some of the false detections that can be produced 
from pipelined image reduction. It should be noted that in 
general for future surveys our transient detection algorithms 
should be capable of recognising (and flagging) common er- 
rors associated with interferometry imaging, hence reducing 
the number of false detections produced. For example, qual- 
ity control measurements should be included in the pipeline 
that assess and extract measurements from the observation 
to remove images from further transient searching. 

Even in a very simplistic implementation, these could 
include, measuring the flux of the calibrator source(s) and 
removing images where the flux had deviated away from the 
correct flux, or ignoring any image where more than > 100 
sources (or any number more than expected) are extracted. 
In this survey, a number of candidate transient sources were 
found to lie on the dirty beam. Source extraction could be 
performed on the dirty beam and checked for associations in 
the CLEANed image. Both false and real transient sources 
can be expected to lie on the dirty beam, however, this in- 
formation could be used to lower the significance of a given 
detection in further analysis. 

If false detections find their way into the database, they 
should be systematically removed to avoid chance source as- 
sociations in future observations. More complex, and com- 



putationally intensive interrogations of the database should 
also be performed to find lower significance detections. 
Greater consideration should also be given to the automated 
flagging algorithms. For this survey flagged data was very 
minimal. For LOFAR, however, large amounts of data may 
be removed, due to the nature of low frequency RFI. Ex- 
amining the flagged visibilities, or even imaging them, could 
yield transient detections. 

This survey did not detect any radio transients and we 
have placed a constraint on the snapshot rate of radio tran- 
sients. We have compared this constraint with results from 
other surveys, and although we did not detect any transients 
it is clear that large volumes of parameter space still remain 
unexplored. As new surveys push into the sub-mjy regime 
with various cadences and at different frequencies, defini- 
tive transient source populations should become apparent. 
Therefore with the next generation of radio telescopes such 
as LOFAR becoming available soon, the problem of inade- 
quate sampling of rare transient phenomena will be allevi- 
ated. Due also to the triggering of multi- wavelength followup 
the classification and interpretation of these events will come 
of age. 

It seems clear that a large population of bright (> 
mjy), GHz, frequent events does not exist. Therefore the 
sub-mjy GHz regime is clearly an important part of pa- 
rameter space to probe for radio transients. The high dy- 
namic range, and wide-field capabilities, of GHz instru- 
ments such as APERTIF (APERture Tiles In Focus, wide 
field upgrade to the W esterbork Synthesis Radio Telescope; 
IVerheiien et al.l 120081 ) and ASKAP make them attractive, 
potentially high yield, discovery instruments. In summary, 
a large population of faint (probably distant) transients, as 
well as very rare bursts remains a strong possibility. Going 
both deeper and wider with the next generation of radio 
facilities will allow us to test these possibilities. 
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